Path-integral Monte Carlo Simulations without the Sign Problem: Multilevel 

Blocking Approach for Effective Actions 
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The multilevel blocking algorithm recently proposed as a 
possible solution to the sign problem in path-integral Monte 
Carlo simulations has been extended to systems with long- 
ranged interactions along the Trotter direction. As an appli- 
cation, new results for the real-time quantum dynamics of the 
spin-boson model are presented. 

PACS numbers: 02.70.Lq, 05.30.-d, 05.40.+j 



I. INTRODUCTION 

Path-integral Monte Carlo (PIMC) simulations are 
useful for extracting exact results on many-body quan- 
tum systems [[[]. In principle, PIMC methods can be 
used to study both equilibrium as well as dynamical 
problems. But in the cases of fermions and real-time 
dynamics, PIMC suffers from the notorious "sign prob- 
lem" which renders such simulations unstable. This sign 
problem manifests itself as an exponential decay of the 
signal-to-noise ratio for large systems or long real times 
[§-f|] . Its origin is at the heart of quantum mechanics it- 
self, namely the interference of different quantum paths 
contributing to the path integral might be destructive 
due to exchange effects or due to the oscillatory nature 
of the real-time evolution operator. Besides approximate 
treatments Q the sign problem has remained unsolved. 

Very recently, a new strategy has been proposed as a 
possible approach to a complete solution of the sign prob- 
lem. This so-called multi-level blocking (MLB) algorithm 
H|J is a systematic implementation of the simple block- 
ing idea — by sampling "blocks" instead of single paths, 
one can always reduce the sign problem 0. Defining a 
suitable hierarchy of blocks by grouping them into dif- 
ferent "levels" , crucial information about the phase can- 
cellations among different quantum paths can then be 
recursively transferred from the bottom to the top level. 
Given sufficient computer memory, such an approach was 
shown to be able to eliminate the sign problem in a sta- 
ble and exact manner ||. But to date, the MLB algo- 
rithm has only been formulated to solve the sign problem 
in PIMC simulations with nearest-neighbor interactions 
along the Trotter direction. This situation is encountered 
under a direct Trotter-Suzuki breakup of the short-time 
propagator. 

In this paper, we report an extension of the MLB ap- 



proach to the case of effective actions that may include 
arbitrarily long-ranged interactions. Such effective ac- 
tions that are non-local in Trotter time may arise from 
degrees of freedoms having been traced out, e.g., a har- 
monic heat bath Q , or through a Hubbard-Stratonovich 
transformation, e.g., in auxiliary- field MC simulations of 
lattice fermions [|3|. Remarkably, because such effective 
actions capture much of the physics, e.g., symmetries or 
the dissipative influence of the traced-out degrees of free- 
dom, the corresponding path integral very often exhibits 
a significantly reduced "intrinsic" sign problem compared 
to the original (time-local) formulation. The present gen- 
eralization of the MLB algorithm was developed to take 
advantage of this fact. We note that in a PIMC simu- 
lation with only nearest-neighbor interactions along the 
Trotter direction, the original MLB approach [|| is more 
efficient than the method described below, which there- 
fore should be used only for time-non-local actions. 

To be specific, we focus on the dynamical sign problem 
arising in real-time PIMC computations here. The mod- 
ifications required to implement the method for fcrmion 
simulations are then straightforward. The structure of 
this paper is as follows. In Sec. |lj the general strategy 
to deal with long-ranged interactions in a MLB scheme 
is outlined. A detailed exposition of the computational 
method can be found in Sec. 
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We have studied the 
real-time dynamics of the celebrated spin-boson system 
H using this approach. Details about this application, 
performance issues related to the sign problem, and nu- 
merical results are presented in Sec. [V. Finally, Sec. |v| 
offers some conclusions. 



II. GENERAL CONSIDERATIONS 

We consider a discretized path integral along a certain 
contour in the complex-time plane. In a typical real- 
time calculation, there is a forward branch from t = 
to t = t* , where t* is the maximum time studied in the 
simulation, followed by a branch going back to the ori- 
gin, and then by an imaginary-time branch from t — 
to t = —iTi(3. We focus on a "factorized" initial prepara- 
tion where the relevant degrees of freedom, denoted by 
r(t), are held fixed for t < That implies that 

the imaginary-time dynamics must be frozen at the cor- 
responding value, and we only need to sample on the two 
real-time branches. Note that such a nonequilibrium cal- 
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culation cannot proceed in a standard way by first doing 
an imaginary-time QMC simulation followed by analytic 
continuation of the numerical data jjj. The quantum 
numbers r(t) at a given time may be discrete or contin- 
uous variables. 

Using time slices of length t*/P, we combine for- 
ward [T"(£ m )] and backward [r'(£ m )] path configurations 
at time t m = mt*/P into the configuration s m , where 
rn = 1,...,P. The configuration at t = is held 
fixed, and for t = t* we must be in a diagonal state, 
r(t*) — r'(t*). For an efficient application of the cur- 
rent method, it is essential to combine several neighbor- 
ing slices m into new "blocks" . For instance, think of 
m = 1,...,5 as a new "slice" £ — 1, m = 6, . . . , 10 
as another slice £ = 2, and so on. Combining q ele- 
mentary slices into a block sg, instead of the original P 
slices we have L = P/q blocks, where L is the number 
of MLB "levels". In actual applications, there is consid- 
erable freedom in how these blocks are denned, e.g. if 
there is hardly any intrinsic sign problem, or if there are 
only few variables in r, one may choose larger values of q. 
Additional flexibility can be gained by choosing different 
q for different blocks. 

Say we are interested in sampling the configurations 
Si on the top level £ — L according to the appropriate 
matrix elements of the (reduced) density matrix, 



p(sl) = Z~ 



E 



exp{-5[si,... > a i ]}, (2.1) 



where S is the effective action under study and Z is a 
normalization constant so that 



Erf 



1 



(2.2) 



Due to the time-non-locality of this action, there will be 



interactions among all blocks sg. The sum in Eq. (2.1) 
denotes either an integration over continuous degrees of 
freedom or a discrete sum. In the case of interest here, 
the effective action is complex-valued and e~ s /\e~ s \ rep- 
resents an oscillatory phase factor (±1 for the fermion 
sign problem). The "naive approach" to the sign prob- 
lem is to sample configurations using the positive definite 
weight function 



7>~|cxp{-S}| , 



(2.3) 



and to include the oscillatory phase in the accumulation 
procedure. Precisely this leads to the exponentially fast 
decay of the signal-to- noise ratio with t* . 

The proposed MLB simulation scheme starts by sam- 
pling on the finest level £ = 1, so only variables in the first 
block corresponding to m = 1, . . . , q are updated. During 
this procedure, interference among different paths will 
take place. Since only relatively few degrees of freedom 
are sampled, however, the resulting interference informa- 
tion can be quantified in a controlled way by employing 



so-called "level-f bonds" (here £ = 1). As long as q is 
chosen sufficiently small, the interference cannot lead to 
numerical instabilities, and the sign cancellations occur- 
ing while sampling on level £ = 1 can thus be synthesized 
and transferred to the level £ = 2, where the sampling 
is carried out next. Here the procedure is repeated, and 
by proceeding recursively up to the top level £ — L, this 
strategy can eliminate the sign problem. The main bot- 
tleneck of the method comes from the immense memory 
requirements, since one needs to store and update the 
level- ^ bonds on all levels during the Monte Carlo sam- 
pling (see below for details). To summarize, the main 
idea of our approach is to subdivide the allowed inter- 
ferences among the quantum paths into small subunits 
(blocks) such that no sign problem occurs when (stochas- 
tically) summing over the paths within each subunit. The 
basic observation underlying our method is therefore al- 
most trivial: The sign problem does not occur in a suffi- 
ciently small system. The nontrivial computational task 
then consists of bringing together the interference sig- 
nals from different blocks, which is done by recursively 
forming blocks on subsequent higher levels. 

Instead of the "circular" structure of the time con- 
tour inherent in the trace operation, it is actually more 
helpful to view the problem as a linear chain, where the 
proposed MLB scheme proceeds from left to right. In 
the case of local actions with only nearest-neighbor inter- 
actions along Trotter time, a different recursion scheme 
was implemented in Refs. |^|| which is close in spirit to 
the usual block-spin transformations used in renormal- 
ization group treatments of spin chains. For both MLB 
implementations, however, the underlying blocking idea 
is identical, and the non-locality of the effective action 
studied here only requires one to abandon block-spin-like 
transformations in favor of the "moving-along-the-chain" 
picture. 

Below we assume that one can decompose the effective 
action according to 



S[Si 



,s L ] =^2Wt[st,...,s L ] 



(2.4) 



All dependence on a configuration sg is then contained 
in the "partial actions" W\ with A < £. One could, of 
course, put all Wi>i = 0, but the approach becomes more 
powerful if a nontrivial decomposition is possible. 



III. MULTILEVEL BLOCKING APPROACH 

In the following, we describe in detail how the MLB 
algorithm for effective actions is implemented in prac- 
tice. The MC sampling starts on the finest level £ = 1, 
where only the configuration sg = i containing the elemen- 
tary slices to = 1, . . . , q will be updated with all S£>i re- 
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maining fixed at their initial values Using the weight 
function 

7>o[«i] - \exp{-W 1 [s 1 ,s 2 ,... > s° L \}\ , 

(i) 

we generate K samples s\, where i — 1,...,K, and 
store them for later use. To effectively solve the sign 
problem and to avoid a bias in the algorithm, the sample 
number K should be chosen large enough, see below and 
Ref . B . For K = 1, the algorithm simply reproduces the 
naive approach. 

The stored samples are now employed to generate in- 
formation about the sign cancellations. All knowledge 
about the interference that occured at this level is en- 
capsulated in the quantity 



B x = 



cxp{-Wi[si, . . . , S L }} 
\exp{-Wi[s u s 2 ) ,...,s L ]}\ / Po[sij 

exp{-W / i[s^ ) , s 2 , . ■ . , s L ]} 



(3.1) 



si 
A" 



U \exp{-W 1 [ S ?\s 2 ,...,sl}}\ 
= Bi[s 2 ,...,s L ] , 

which we call "level- 1 bond" in analogy to Ref. ||, with 
the normalization constant C'o — E S i "Po[ s i\- The third 
line follows by noting that the were generated ac- 
cording to the weight Vq. This equality requires that K 
is sufficiently large and that q is sufficiently small in or- 
der to provide a good statistical estimate of the level- 1 
bond. 



Combining the second expression in Eq. (3.1) with 
Eq. ( |2.lD , we rewrite the density matrix in the follow- 
ing way: 

p(s L ) = Z- 1 expl-J2wAc B 1 (3.2) 

S 2 ,---,SL-1 I i>l ) 



z- 1 ]T PoSiJJ. 



«1 «L-1 



1>1 



When comparing Eq. (3.2) with Eq. (2.1), we see that 



the entire sign problem has now formally been trans- 
ferred to levels I > 1, since oscillatory phase factors only 
arise when sampling on these higher levels. Note that 
B\ = Bx[s2-, ■ ■ ■ , sl] introduces couplings among all lev- 
els £ > 1, in addition to the ones already contained in 
the effective action S. 

We n ow p roceed to the next level I = 2 and, according 
to Eq. (3.2), update configurations for m = q + 1, . . . , 2q 
using the weight 

Vi[a 2 ] = \B 1 [ S2 ,s°,...,sl}eM-W2[s2,sl...,s° L }}\ . 

(3.3) 



Under the move s 2 — > s' 2 , we should then resample and 
update the level-1 bonds, B\ — ► B' x , Exploiting the fact 
that the stored K samples are correctly distributed 
for the original configuration s 2 , the updated bond can 
be computed according to 



K 

X E 



exp{-W 1 [s^,s' 2 ,...,s L }} 



~{ \exp{-Wi[s?>,8 



21 ■ 



(3.4) 



Again, to obtain an accurate estimate for B[ , the number 
K should be sufficiently large. In the end, sampling under 
the weight V\ implies that the probability for accepting 
the move s 2 — > s' 2 under the Metropolis algorithm is 



cxp{-W 1 [ B ^,^,^,...]} 
,a°,...]}\ 



E, 



C xp{~W 1 [s[ l) ,s 2 ,s^...}} 



| cxp{ — Wi [s 



.§,...]}! 



exp{-W 2 [s 2 ,sl...}} 



exp{-W / 2 [s2, s°, • 



Using this method, we generate K samples s 2 \ store 
them, and compute the level-2 bonds, 

o = / B 1 [s 2 ,s 3 ,...]exp{-W 2 [s2,s 3 ,...]} \ 
2 \\B 1 [s 2 ,sl...}exp{~W 2 [s 2 ,sl...}}\/ Vi[s2] [ -> 



]} 



(3.5) 



Cf^Sitaa,...] exp{-W 2 [s 2 ,...}} 



K 



^E 



s 3 ,...]exp{-W 2 [s 2 z \s 3 ,...]} 



=T |B 1 [4 <) )S o,...]exp{-W 2 [^ )S ; 



2 >°3i' 



]}| 



B 2 [s 3 , . ..,S L ] , 



with Ci = J2s 2 ^ 3 i[ s 2]- Following our above strategy, we 
then rew rite the reduced density m atrix by combining 
Eq. (3^) and the second line of Eq. (3J3). This yields 



Z~ 



= Z~ 



E 



exp 



l>2 



W t } C dB 2 (3.7) 



Si,...,S_l_i 



e>2 



Clearly, the sign problem has been transferred one block 
further to the right along the chain. Note that the nor- 
malization constants Co, C\, . . . depend only on the initial 
configuration s° so that their precise values need not be 
known. 

This procedure is now iterated in a recursive manner. 
Sampling on level £ using the weight function 

Ve-i[s e ] = |S £ _i[s£, s° +1 , . . .} exp{-W e [s e , s° e+1 , . . .]}| 

(3.8) 

requires the recursive update of all bonds B\ with A < 
I. Starting with B\ — > B[ and putting Bq = 1, this 
recursive update is done according to 
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B[ = K~ 



E 



(3.9) 

exp{-^[4 l) , SA+1 ,...]} 



i \Bx-x[sf,s» 



'A+l> 



exp{-WA[4 l) , S ° +1 , 



•]}| 



where the primed bonds or partial actions depend on s' e 
and the unprimed ones on s°. Iterating this to get the 



updated bonds Bg-2 for all s^_ v the test move — > s', 

(3.10) 



is then accepted or rejected according to the probability 
Bt^x[s' t , s% +1 , . . .} cxp{~W e [s' e , s° e+1 , . . .]} 



Bt-i[s£, s% +1 , . . .} exp{-We[s e , s% +v . . .]} 



(i) 

On this level, we again generate K samples s\ , store 
them and compute the level-£ bonds according to 



B e [s i+1 , ...}= K 1 



(3.11) 



K 
i=l 



Bt^i[s^\s t+ i, . . .] exp{-Wi[af , s e+ i, . . .]} 



, r„W s o 



|exp{-Wi[a^,«S +1 ,...]} 



This process is iterated up to the top level, where the 
observables of interest may be computed. 

Since the sampling of Bi requires the resampling of all 
lower-level bonds, the memory and CPU requirements of 
the algorithm laid out here are quite large. For A < £—1, 
one needs to update B\ — > B' x for all sf) with A < £' < £, 
which implies a tremendous amount of computer mem- 
ory and CPU time, scaling approximately ~ K L at the 
top level. Fortunately, an enormous simplification can 
often be achieved by exploiting the fact that the inter- 
actions among distant slices are usually weaker than be- 
tween near-by slices. For instance, when updating level 
£ = 3, the correlations with the configurations may 
be very weak, and instead of summing over all K samples 

(i) 

s\ in the update of the bonds -Ba<£, we may select only 
a small subset. When invoking this argument, one should 
be careful to also check that the additional interactions 
coming from the level-A bonds with A < £ are sufficiently 
short-ranged. From the definition of these bonds, this is 
to be expected though. 

Remarkably, this algorithm can significantly relieve the 
severity of the sign problem. Let us first give a simple 
qualitative argument supporting this statement for the 
original MLB method of Ref. ||, where P = 2 L with 
L denoting the number of levels. If one needs K sam- 
ples for each slice on a given level in order to have sat- 
isfactory statistics despite of the sign problem, the total 
number of paths needed in the naive approach depends 
exponentially on P, namely ~ K p . This is precisely the 
well-known exponential severity of the sign problem un- 
der the naive approach. However, with MLB the work 
on the last level [which is the only one affected by a 
sign problem provided K was chosen sufficiently large] 
is only ~ K L . So in MLB, the work needed to sample 
the K p paths with satisfactory statistical accuracy grows 



~ K l °zi p = P lo «2 K i i.e., only algebraically with P. Pro- 
vided the interactions along the Trotter time decay suffi- 
ciently fast, a similar qualitative argument can be given 
for the generalized MLB algorithm proposed here. For 
the application described below, we have indeed found 
only algebraic dependences of the required CPU times 
and memory resources with the maximum real time t* , 
instead of exponential ones as encountered in the naive 
approach. Further details of the simulation procedure 
are provided in the next section. 



IV. APPLICATION: SPIN-BOSON DYNAMICS 

To demonstrate this MLB algorithm for path integral 
simulations with long-range interactions in the Trotter 
direction, we study the real-time dynamics of the spin- 
boson model, 



H 



-(HA/2) tr x 

,2 



(he/2) a z 



(4.1) 



E 



Pn 



2m r 



\m a uj 2 a 



This archetypical model has a number of important ap- 
plications, e.g., the Kondo problem, interstitial tunnel- 
ing in solids quantum computing pXfl , and electron 
transfer reactions pj] , to mention only a few. The bare 
two-level system (TLS) has a tunneling matrix clement 
A and the asymmetry (bias) e between the two localized 
energy levels (a x and a z are Pauli matrices). Dissipation 
is introduced via a linear heat bath, i.e., an arbitrary col- 
lection of harmonic oscillators {x a } bilinear ly coupled to 
a z . Concerning the TLS dynamics, all information about 
the coupling to the bath is contained in the spectral den- 
sity J(oj) = (ir /2)J2 a ( c a/ m a UJ a) S(u> — w a ), which has a 
quasi-continuous form in typical condensed-phase appli- 
cations. J(u>) dictates the form of the (twice-integrated) 
bath correlation function (/3 = l/fc^T), 



Q(t) = 



duj JIuj) cosh[w?i/?/2] - cosh[w(fi,/3/2 - it)} 



smh[uhl3/2] 



(4.2) 



For the calculations here, we assume an ohmic spectral 
density of the form J(jjj) — 2ttTicvjJ exp(— u>/u> c ), for which 
Q(t) can be found in closed form [Q. Here lu c is a cutoff 
frequency, and the damping strength is measured by the 
dimensionless Kondo parameter a. In the scaling limit 
A -C w c , and assuming a < 1, all dependence on uj c 
enters via a renormalized tunnel splitting S 



ieff 



[cos(7ra)r(l - 2a)] 1/2(1 - Q) (A/w c ) Q/(1 - 



(4.3) 



and powerful analytical |8|,|T2| anc ^ alternative numerical 
methods [ |l3|Jl4| ] are available for computing the nonequi- 
librium dynamics. 
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At this point some remarks are in order. Basically 
all other published numerical methods except real-time 
PIMC can deal only with equilibrium quantities, see, 
e.g., Refs. [15,16], or explicitly introduce approxima- 
tions |EMMl7,lS]. Regarding the latter class, mostly 
Markovian-type approximations concerning the time- 
range of the interactions introduced by the influence func- 
tional have been implemented. Our approach is computa- 
tionally more expensive than other methods [fL3|-p"8| , but 
at the same time it is unique in yielding numerically exact 
results for the nonequilibrium spin-boson dynamics for 
arbitrary bath spectral densities. It is particularly valu- 
able away from the scaling regime where important appli- 
cations, e.g., coherent (nonequilibrium) electron transfer 
reactions in the adiabatic regime, are found but basically 
all other methods fail to yield exact results. Finally we 
briefly compare the present approach to our previously 
published PIMC method 0. For not exceedingly small 
a, it turns out that the latter method is just equivalent 
to the K = 1 limit of the present method. From Table I 
and the discussion below, it is thus apparent that MLB 
is significantly more powerful in allowing for a study of 
much longer real times than previously. 

We study the quantity P(t) = {<y z {t)) under the 
nonequilibrium initial preparation a z (t < 0) = +1. 
P{t) gives the time-dependent difference of the quantum- 
mechanical occupation probabilities of the left and right 
states, with the particle initially confined to the left state. 
To obtain P(t) numerically, we take the discretized path- 
integral representation of Ref. |?J and trace out the bath 
to get a long-ranged effective action, the "influence func- 
tional" . In discretized form the TLS path is represented 
by spins <Ji,<j[ = ±1 on the forward- and backward-paths, 
respectively. The total action S consists of three terms. 
First, there is the "free" action So determined by the 
bare TLS propagator Uq, 

p-i 

exp(-So) = J] U (a i+1 ,ai;t*/P) U (a' i+1 , a' t ; -t*/P) . 

(4.4) 

The second is the influence functional, Si = Sj + Sj , 
which contains the long-ranged interaction among the 
spins, 

) = E - °j){ L J-m - <4) (4-5) 

j>m 

+ iL j-rn ( CT m + O} ) 

where Lj = L'j + iLj is given by 
Lj = [Q((j + l)t*/P) + Q((j - l)t*/P) - 2Q(jt*/P)}/4 

(4.6) 

for j > 0, and Lq = Q(t* / P)/4. In the scaling regime 
at T = 0, this effective action has interactions ~ a/t 2 



between the spins ("inverse-square Ising model"). The 
contribution 



Sf } = i(f/P) E 7K /P)(a m a'J 



(4.7) 



gives the interaction with the imaginary-time branch 
[where a z = +1], where the damping kernel 



7(*) = 



■kTi 



OLD COS(ojt) 



(4.8) 



For clarity, we focus on the most difficult case of an un- 
biased two-state system at zero temperature, e = T = 0. 
To ensure that the Trotter error is negligibly small, we 
have systematically increased P for fixed t* until conver- 
gence was reached. Typical CPU time requirements per 
10 4 MC samples are 4 hours for P = 26, L = 2, K = 1000, 
or 6 hours for P — 40, L = 3, K — 600, where the simula- 
tions were carried out on SGI Octane workstations. The 
memory requirements for these two cases are 60 Mbyte 
and 160 Mbyte, respectively. Data were collected from 
several 10 5 samples. 

For a = 0, the bare TLS dynamics P(t) — cos(At) is 
accurately reproduced. As mentioned before, the perfor- 
mance is slightly inferior to the original MLB approach 
Q which is now applicable due to the absence of the in- 
fluence functional and the associated long-ranged inter- 
actions. Turning to the situation where a bath is present, 
we first study the case a = 1/2 and w c /A = 6. The exact 
a = 1/2 result ||, P(t) = cxp(— A c gt), valid in the scal- 
ing regime w c /A 3> 1, was accurately reproduced, indi- 
cating that the scaling regime is reached already for mod- 
erately large u c /A. Typical parameters used in the MLB 
simulations and the respective average sign are listed in 
Table |[ The first line in Table | corresponds to the naive 
approach. For a = 1/2, it turns out that our previous 
PIMC scheme [[Fj yields a comparable performance to the 
K = 1 version of this MLB method. It is then clear from 
Table Q that the average sign and hence the signal-to- 
noise ratio can be dramatically improved thus allowing 
for a study of significantly longer timescales t* than be- 
fore. For a fixed number of levels L, the average sign 
grows by increasing the parameter K. Alternatively, for 
fixed K, the average sign increases with L. Evidently, 
the latter procedure is more efficient in curing the sign 
problem, but at the same time computationally more ex- 
pensive. In practice, it is then necessary to find a suitable 
compromise. 

Figure shows scaling curves for P(t) at a = 1/4 for 
u>c/A = 6 and ui c /A = 1. According to the a = 1/2 
results, w c /A — - 6 is expected to be within the scaling 
regime. This is confirmed by a comparison to the non- 
interacting blip approximation (NIB A) Q. The minor 
deviations of the NIBA curve from the exact result are 
in accordance with Refs. j712| for a < 1/2. However, 
for oj c /A = 1, scaling concepts (and also NIBA) are ex- 
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pected to fail even qualitatively. Clearly, the MLB re- 
sults show that away from the scaling region, quantum 
coherence is able to persist for much longer, and both 
frequency and decay rate of the oscillations differ signif- 
icantly from the predictions of NIBA. In electron trans- 
fer reactions in the adiabatic-to-nonadiabatic crossover 
regime, such coherence effects can then strongly influence 
the low-temperature dynamics. One obvious and impor- 
tant consequence of these coherence effects is the break- 
down of a rate description, implying that theories based 
on an imaginary-time formalism might not be appropri- 
ate in this regime. A detailed analysis of this crossover 
regime using MLB is currently in progress. 

V. CONCLUSIONS 

In this paper, we have extended the multilevel block- 
ing (MLB) approach of Refs. (|]](| to path-integral 
Monte Carlo simulations with long-ranged effective ac- 
tions along the Trotter direction. For clarity, we have 
focussed on real-time simulations here, but believe that 
a similar approach can also be helpful in many-fermion 
computations, e.g., in auxiliary-field fermion simulations 
of lattice fermions. The practical usefulness of the ap- 
proach was demonstrated by computing the nonequi- 
librium real-time dynamics of the dissipative two-state 
system. Here the effective action (influence functional) 
arises by integrating out the linear heat bath. For a heat 
bath of the ohmic type, at T = the corresponding in- 
teractions among different time slices decay only with a 
slow inverse-square power law. 

In the present implementation of MLB, the basic block- 
ing idea operates on multiple time scales by carrying out 
a subsequent sampling at longer and longer times. Dur- 
ing this procedure, the interference information collected 
at shorter times is taken fully into account without in- 
voking any approximation. Under such an approach, at 
the expense of large memory requirements, the severity 
of the sign problem can be significantly relieved. The 
proposed approach allows to study time scales not acces- 
sible to previous real-time path-integral simulations for 
the spin-boson system. 
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FIG. 1. Scaling curves for P(t) for a = 1/4 with lj c /A = 6 
(closed diamonds) and cu c /A = 1 (open circles). The solid 
curve is the NIBA prediction. The approach of Ref. [jjj be- 
comes unstable for A e a£ > 4 in both cases. Statistical errors 
are of the order of the symbol sizes. 

TABLE I. MLB performance for a = 1/2, to c /A = 6, 



At* = 10, P = 40, and several L. qe denotes the number of 
slices for I = 1, . . . L. 

K_ L (sgn) 

1 1 40 0.03 

200 2 30 - 10 0.14 

800 2 30 - 10 0.20 

200 3 22 - 12 - 6 0.39 

600 3 22 - 12 - 6 0.45 
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